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ABSTRACT 


— This  study,  assuming  stationary  Gaussian  turbulence 
model,  investigates  the  effect  on  the  crossing  frequency  of 
different  spectral  functions:  von  Karman,  Kaimal  and 
Teunissen  for  the  x  and  z  directions  with  five  filters: 
ideal  band-pass,  ideal  low-pass,  "quadra tic-type",  "sine- 
type”  and  the  Hanning.  The  filters  have  a  much  greater 
effect  than  the  spectral  functions.  The  estimated  crossing 
frequency  variation  is  as  much  as  50  percent  jmong  the 
quadratic-type,  low-pass  and  sine-type  filters.  The  Hanning 
filtering  predicts  crossing,  rates  up  to  thirty  times,  and 
the  ideal  one  decade  wide  band-pass  filtering  predicts  rates 
of  eight  times,  higher  than  the  ideal  low-pass  filtering. 
The  variation,  between  the  x  and  z  direction,  is  less  than 
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I.  INTRODUCTION 


Most  flows  occurring  in  nature  and  in  engineering  appli¬ 
cations  are  turbulent.  In  the  earth's  boundary  layer, 
differential  heating  of  the  atmosphere  produces  pressure 
gradients,  which  are  subsequently  modified  by  the  rotation 
of  the  earth,  causing  a  complex  velocity  field.  In  this 
boundary  layer  the  wind  speed  decreases  as  the  surface  is 
approached  due  to  both  the  frictional  drag  of  the  surface 
and  the  drag  of  all  bodies  protruding  into  the  air  flow. 
These  retarding  forces  are  transmitted  through  the  layer  by 
shear  forces  and  by  the  exchange  of  momentum  due  to  the 
vertical  movement  of  the  air.  The  process  of  momentum 
exchange  between  layers  is  the  mechanism  leads  to  the  gener¬ 
ation  and  decay  of  eddies  which  are  termed  turbulence.  The 
resulting  mixing  of  the  air  produces,  along  all  three 
orthogonal  axes,  fluctuations  in  wind  speed,  commonly  called 
gusts,  which  vary  in  size  in  both  time  and  space. 

The  effects  of  this  atmospheric  turbulence  has  been  of 
continuing  concern  to  the  aircraft  or  structure  designer. 
Typical  turbulence  related  problems  are:  the  effects  of 
turbulence  on  the  fatigue  life  of  the  structure;  the 
performance  of  control  systems  in  turbulence  and  the  deter¬ 
mination  of  ultimate  structural  strength  required  to  accept 
unsteady  loads  induced  by  turbulence.  In  an  attempt  to  solve 
these  problems,  a  number  of  statistical  models  of  turbu¬ 
lence,  which  endeavour  to  describe  the  turbulence  in  terms 
of  as  few  parameters  as  possible,  have  been  proposed  since 
early  in  this  century.  To  analyze  these  problems,  one  needs 
to  know  the  crossing  frequency,  which  is  the  frequency  with 
which  a  random  function (  here  the  filtered  wind  speed) 
crosses  a  prescribed  value  defined  by  a  certain  averaging 
time,  say  the  hourly  mean. 


The  crossing  frequency  depends  largely  on  whether  or  not 
the  random  function  is  filtered,  the  type  of  filter  and  to  a 
lesser  extent, the  spectral  density  function.  The  effects  of 
type  of  filter  used  and  spectral  density  expression  will  be 
discussed  in  detail  later.  The  spectral  density  functions 
have  been  proposed  by  Dryden  [Ref.  1],  von  Karman  [Ref.  2], 
Kaimal  [Ref.  3]  and  Teunissen  [Ref.  4],  The  Dryden  expres¬ 
sion  makes  mathematical  analysis  simple  but  is  much  less 
accurate  than  the  others.  The  von  Karman  spectral  function 
was  recommended  by  the  Engineering  Sciences  Data  Unit 
(E.S.B.O.)  [Ref.  5]  and  this  expressions  was  originally 
postulated  for  isotropic  turbulence.  The  effect  of  departure 
from  isotropic  turbulence  near  the  ground  is  allowed  for  by 
the  variation  with  height  and  surface  roughness  of  the 
appropriate  variance  and  length  scale  parameter  since 
they  typify  the  intensity  and  size  of  eddies  constituting 
turbulence.  The  Kaimal  expression,  which  is  frequently 
used,  was  obtained  recently  for  the  surface  layer  over  flat, 
uniform,  relatively  featureless  terrain  in  Kansas.  The 
Teunissen  model,  which  is  a  modified  Kaimal  one,  was 
obtained  for  the  generally  rougher  gross  features  of  the 
upstream  terrain. 

In  the  following  discussion,  we  will  review  turbulence 
models  and  spectral  density  functions.  The  spectral  density 
S(n),  in  non-dimensional  form,  is  a  function  of  the  dimen¬ 
sionless  parameter  nL^/O,  where  n  is  frequency  and  Lj  is  the 
length  scale  of  the  turbulence  and  U  is  the  mean  hourly  wind 
speed.  E.S.D.U.  [Ref.  5]  gives  empirical  expressions  for 
the  length  scale  obtained  by  analysing  its  collection  of 
world-wide  turbulence  data.  Finally,  for  given  values  of 
nL^/0,  we  will  estimate  the  crossing  frequencies  with  von 
Karman,  Kaimal  and  Teunissen  spectral  expressions. 
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II.  TURBULENCE  MODELS 


A.  DISCRETE  GOSTS  MODEL 

Powell  and  Connell  £Ref.  6]  defined  gusts  as  constitu¬ 
ting  any  series  of  discrete  velocity-time  events  that  can  be 
defined  from  a  turbulence  time  series  according  to  some 
extrinsic  criterion.  There  are  two  fundamentally  different 
treatments  of  gust  time;  one  is  that  the  gust  time  is  arbi¬ 
trarily  fixed  and  the  other  is  that  the  gust  time  varies. 
Examples  of  this  latter  type  of  gust  definition,  shown  in 
Figure  2.1,  represent  a  sample  of  wind  fluctuation  in  the 
form  of  a  wind-component  time  series. 

The  GUSTO  and  GUST1  models  define  gust  events  in  terms 
of  gust  amplitudes  and  characteristic  times.  Definitions  of 
amplitudes  and  times  differ  for  the  two  models.  The  GUSTO 
model  can  be  completely  specified  in  terms  of  a  positive 
peak  amplitude  +  A0  ,  and  the  time  interval  *T0  between  zero 
crossings  on  either  side.  It  can  also  be  specified  in  terms 
of  a  negative  peak  amplitude  -A0  and  the  corresponding  time 
interval  ~T0  .  Both  definitions  are  expected  to  yield 
similiar  statistics.  Powell  and  Connell  combine  both  posi- 
-e  and  negative  values  into  one  set  for  their  analysis. 

The  GUST1  model  can  be  specified  in  terms  of  a  positive 
amplitude  +At  which  is  the  peak-to-peak  amplitude  between 
adjacent  minima  and  maxima,  and  a  time  interval  +T,  between 
the  minima  and  the  maxima;  the  positive  sign  indicates  a 
positive  rate  of  change  in  the  variable.  A  comparable  defi¬ 
nition  can  be  made  in  terms  of  amplitude  -Aj  and  time  -Tt  , 
where  the  negative  sign  indicates  a  negative  rate  of  change 
in  the  variable. 
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One  can  find  the  impulse  function  of  Equation  (4. 9)  by 
inverse  Fourier  transform/  i.e., 

h4(t)  =7rjjnc  +n2  ),/2ej21cntdn 

=  I (n2  +n2  )'^ccs2rtntdn  + 

-^-j(n*  +  n2  )l/a  ( jsin2rcnt)  dn  (4-  10) 

The  second  term  of  above  equation  is  zero,  since,  it  is  an 
odd  function.  Thus  Equation (4 .  10)  becomes 

-  <70 

h4(t)  =^-  J  (n*  +n2)^  cos2i intdn 
Jo 

or  h^(t)  =  -£-K0(nct)  (4.11) 

where  K0 is  the  modified  Bessel  function. 

The  functions  h^  (t)  and  F/j.  (n)  are  shown  in  Figure  4.1 
and  Figure  4.2. 

The  cut-off  frequencies  are  defined  formally  as  the 
frequencies  at  which  the  filters  reach  the  first  zero.  But, 
when  a  real  filter  is  used,  we  must  find  some  equivalent 
cut-off  frequency,  which  is  defined  as  the  cut-off  frequency 
of  an  equivalent  ideal  filter  whose  area  is  equal  to  that  of 
the  real  filter.  [Bef.  11]  Greenway  [Ref.  12]  and  Powell 
and  Connell  [Bef.  6]  have  used  the  value  at  the  half- 
power  (3db)  point  of  the  filter  as  the  cut-off  frequency. 
Vinnichenko  et  al.  [Ref.  11]  defined  the  equivalent  cut-off 
frequency  as 


n 


c 


(4.  12) 
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FILTER  FUNCTION 


The  corresponding  transfer  function  F3  (n)  has  the  form 

F3<n)  =  (  sinitTn^2  (4.8) 

\  TCTn  I  [lU- (TTTn)1]2, 

The  advantage  of  F3  (n)  over  F-^n),  is  that  it  has  much 
smaller  side  maxima. 

Another  commonly  used  filter  function  is  the  guadratic 
filter  function  given  by 

(n)  =  1  (4.9) 

n*  ♦  n2. 
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TIME  FUNCTION 


Figure  4.1  Tine  Function  ht  (t)  ,  h2  (t)  ,  h3(t)  and  h+(t)  . 


The  function  h2  (t)  and  F2  (n)  are  shown  in  Figure  4.  1  and 
Figure  4.2  respectively.  The  figure  shows  that  Fz  (n)  is  a 
low-f reguency  filter  with  cut-off  frequency  nc  =1/T.  The 
shape  of  this  filter  shows  that  frequencies  above  nc  are 
transmitted  nonuniformly  and  those  below  nc  are  completely 
suppressed. 

The  impulse  response  function  h3  (t) ,  called  the  smoothed 
moving  average  is: 


b3  (t) 


|  1«-cos{2*/T)  t 


when 


-T/2  <t<  T/2 


(4.7) 


0 


when  1 1 1  >  T/2 


A.  FILTE8  FUNCTIONS 


It  is  clear  that  the  quality  of  the  crossing  frequency 
estimates  obtained  by  means  of  Equation  (4. 2)  will  depend 
strongly  on  the  form  of  the  filter  function  F(n)  and,  hence, 
on  the  selection  of  the  time  function  h(t).  Thus  we  will 
consider  some  examples  of  the  function  h  (t)  and  filter  func¬ 
tion  F  (n)  corresponding  to  them.  Let  the  time  function  be 
given  by 


hL(t)  = 

2  sin  (2xt/T) 

(4.3) 

T  2ltt/T 

If  we  take  Fourier  transform. 

t  he 

Equation  (4.3) 

becomes 

Fj.  (n)  =| 

r  1  when 

l®l 

<  1/T 

(4.4) 

1 

1/2  when 

mi 

=  1/T 

.  0  when 

mi 

>  1/T. 

Graphs  of  the  functions  hj  (t)  and  Ft  (n)  are  shown  in  Figure 
4.1  and  Figure  4.2  respectively  and  they  represent  an  ideal 
band-pass  filter. 

The  simple  impulse  response  function  h2  (t)  is  the  ordi¬ 
nary  moving  average: 

h2  (t)  1/T  when  -T/2  £t<  +T/2  (4.5) 

1 0  when  |  t {  >  T/2 


The  corresponding  filter  function  F2  (n)  (F2  (n)=H2(n)'i  ),  where 
H2(n)  is  the  Fourier  transform  of  the  function  h2(t),  has 
the  form 


F2  (n) 


sinxTn 

ltTn 


(4.6) 
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17.  CROSSING  PEEQOENCIES 


The  instantaneous  horizontal  wind  speed  is  considered  to 
be  composed  of  a  mean  longitudinal  component  0  plus  fluctu¬ 
ating  component  v;  ,  where  the  subscript  i  designates  the  x, 
y,  or  z  component  in  a  coordinate  system  with  x  oriented 
along  the  mean  wind  vector  and  z  upwards.  Then  the  wind 
speed  as  a  function  of  time  is  given  by 


V  (t)  =  0  +  V;  (t) 


(4.1) 


The  average  frequency  of  positive-slope  level  crossing 
Noi ,  with  which  the  wind  speed  fluctuations  exceed  their 
zero  value,  for  a  Gaussian  distribution  of  turbulence,  is 
given  by 


Noi  *  I  ^FfnJS;  (n)dn  /  (  F  (n)  St-  (n)  dn 


(4.2) 


where  nc  is  cut-off  freguency,  i  is  u,  v,or  w  and  F  (n)  is  a 
filter  function. 

The  value  of  N0^*  depends  on  the  mean  wind  speed  U,  Lf  and 
turbulence  intensity  c£/U,  which  is  a  measure  of  the  magni¬ 
tude  of  turbulence  fluctuations  and  defined  as  the  ratio  of 
the  standard  deviation  of  the  instantaneous  fluctuating 
velocity  component  to  the  mean  wind  speed  averaged  over 
arbitrary  time,  usually  hourly,  and  depends  only  on  the 
height  z  and  the  terrain  roughness  z0  ,  in  a  stable 
atmosphere. 


function  of  the  gust  velocity  can  be  found  by  integration. 


i.e. , 


r<? 

P  =  /  p  (i)  di  (3.9) 

-^.00 

In  practice,  atmospheric  turbulence  contains  patches  of 
a  significantly  non-Gaussian  nature  (particularly  in  the 
lower  30m)  when  larger  gusts  and  longer  lulls  may  occur  more 
frequently  then  indicated  by  the  Gaussian  distribution. 
Kaimal  et  al.  £Hef.  10]  observed  that  the  assumption  of  a 
stationary  Gaussian  process  has  validity  only  for  the 
filtered  data,  not  for  the  unfiltered  time  series. 
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where  f=nL^/U  and  u*  is  friction  velocity.  These  expres¬ 
sions  were  derived  from  the  data  which  were  obtained  from  a 
horizontal  array  of  tower-mounted  propeller  anemome¬ 
ters  (z=11m)  during  a  five-hour  period  for  which  the  mean 
wind  direction  was  virtually  perpendicular  to  the  main  span 
of  the  array,  in  the  city  of  Toronto,  Canada. 

As  for  the  reason  for  the  higher  spectral  energies  found 
by  Teunissen  in  the  low-f reguency  region  of  the  spectra  in 
comparison  with  the  Kaimal  values  which  was  measured  in 
Kansas,  this  is  not  entirely  obvious,  apart  from  the  rougher 
terrain  for  Teunissen  experiment.  Teunissen  has  suggested 
that  the  power  spectra  over  a  wide  range  of  terrain  types 
may  possibly  be  represented  by  using  above  expressions  with 
appropriate  values  of  such  ’terrain  scaling’  parameters. 

B.  PBOBABILITY  DENSITY 

In  general,  a  knowledge  of  the  spectral  density  of  a 
random  process  does  not  enable  us  to  determine  its  prob¬ 
ability  density  or  distribution.  However,  a  common  assump¬ 
tion  which  is  reasonable  in  many  applications  is  that 
atmospheric  turbulence  is  a  ’normal’  or  Gaussian  process,  as 
discussed  above,  with  a  probability  density  function  for 
which 


P(i)'*/Tir  “eh1''2'*1]  (3-8) 

where  i=u(t),  v(t)  or  w(t)  are  wind-speed  fluctuations, 
which  have  a  zero  mean.  Note  that  Equation  (3. 8)  depends  only 
on  the  standard  deviation  of  the  gust  velocity.  Thus,  if  the 
probability  distribution  of  <r  is  known,  the  distribution 


Another  commonly  used  spectral  expression  is  the  Kaimal 
one;  [Ref.  3] 

nSj  (n)  _  0.164  (f/f«)  (3.4) 

c^a  1+0.  164{f/f0  )5/a 

where  nS  (n)  =  logarithmic  power  spectral  density 
d"l  -  variance  of  i 

f  =  reduced  freguency 

f„  =  reduced  freguency  at  the  intercept  of  the 
extrapolated  inertial  subrange  slope  with  the  nS4*(n)/ 
O'/ =1  line;  or  in  our  notation,  f*;  =0. 041z/L4*  . 

The  Kaimal  expression  is  acre  common  in  the  meteorolo¬ 
gical  literature  and  was  obtained  as  a  best-fit  to  surface- 
layer  measurements  over  uniform,  flat,  relatively 
featureless  terrain. 

Teunissen  [Ref.  4]  observed  that  von  Karman  model 
appears  to  be  somewhat  better  than  the  Kaimal  model  at  low 
freguencies.  The  Kaimal  model  largely  underestimates  the 
spectral  content  at  these  freguencies  for  all  compo¬ 
nents  (x,  y,  and  2).  On  the  other  hand,  the  Kaimal  model  shape 

appears  to  be  better  than  that  cf  the  von  Karman  model. 

Teunissen  has  proposed  modified  Kaimal  spectral  expres¬ 
sions,  given  by 


nSw  (n)  _ 

1 0  5f 

(3.5) 

(0.44*33f)5/5 

nSv(n)  _ 

1  7£ 

(3.6) 

c 

1 

(0.38*9.5f//3 

nS«/(n)  _ 

2  f 

(3.7) 

C 

1 

0.44*5.3f5/3 
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Spectral  functions  of  atmospheric  turbulence  provide 
information  on  the  frequency  distribution  of  the  kinetic 
energy  of  the  various  fluctuating  velocity  components.  Used 
in  conjunction  with  certain  transfer  functions,  they  provide 
information  about  the  dynamic  loading  on,  and  response  of, 
building  and  aircraft  structures  in  the  atmospheric  wind. 

A  considerable  number  of  measurements  of  the  power 
spectra  of  gust  velocities  in  the  near-neutral  atmosphere, 
or  strong  wind  conditions,  at  varying  heights  and  for 
different  terrains  are  available.  The  von  Karman  spectral 
expressions,  which  were  exclusively  used  by  E. S.D.U. 
[fief.  5],  for  each  velocity  component  (x,  y  and  z  direction) 
are  given  by 


nSu  (n)  _  4nu  (3.1) 

<Tu  (1  +  70.8ni)5/6 


nSy  (n)  ___  4s y  (1^7 55.  2n  y )  (3.2) 

cr$  ( 1+283.  2nv)"/6 


nSw  (n)  _  4nm  (1+755.  2n£)  (3.3) 
crj  (1+28  3.255)*^ 

where  nu=  Lun/U,  nv=  L„n/U,  nw=  L* n/U  and  Lp  is  integral 
length  scale  which  is  a  function  only  of  the  terrain  rough¬ 
ness  and  height  above  ground  for  the  x  and  y  directions;  for 
the  z  direction  it  is  a  function  of  elevation  only,  o'/  is 
the  variance  of  the  wind  velocity  fluctuations  about  the 
hourly  mean.  These  equations  are  in  common  use  in  the  engi¬ 
neering  literature  and  have  the  advantage  that  the  integral 
length  scale  are  treated  as  'free'  scaling  parameters 
which  are  chosen  to  match  the  estimated  scales  for  a  partic¬ 
ular  height  and  terrain  type,  while  maintaining  constant 
spectral  shape. 
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III.  SgECTfi&L  DEHSITY  EROB^BILITX  DENSITY 


A.  SPECTRAL  DENSITY 

In  the  field  of  mechanical  and  electrical  oscillations, 
a  large  number  of  phenomena  are  governed  by  equations  of  the 
form 


■  fe  *  c  *  **  ’  £(t) 

where,  m,  c  and  k  are  constants;  f  (t)  is  a  given  forcing 
function  and  y (t)  represents  the  response  of  the  system  as  a 
function  of  the  time  t. 

If  the  above  time  dependent  force  f  (t)  is  periodic,  it 
can  be  expressed  in  terms  of  a  Fourier  series;  this  series 
represents  the  sum  of  a  large  number  of  forces  of  different 
amplitudes  and  frequencies.  All  elastic  bodies  have  natural 
frequencies,  and  if  this  force  is  applied  to  an  elastic 
system,  and  the  amplitudes  of  the  components  near  a  natural 
frequency  are  not  very  small,  then  this  force  will  drive  the 
body  into  resonance.  The  body  will  then  go  through  a  very 
large  number  of  stress  cycles  in  a  short  period  of  time, 
possibly  leading  to  fatigue  failure. 

To  have  stationary  properties,  a  random  signal  must  be 
assumed  to  continue  over  an  infinite  time,  and  in  such  a 
case  neither  the  real  nor  the  imaginary  part  of  the  Fourier 
transform  converges  to  a  steady  value.  In  this  case,  we  use 
the  spectral  density  function,  which  has  no  convergence 
difficulty  and  which  is  applicable  to  a  whole  class  of 
similarly  generated  functions. 


behave  like  exp(-x)  rather  than  exp(-x  ),  as  predicted  by 
the  Gaussian  distribution  [Ref.  8].  Recently  Reeves,  et  al. 
[Ref.  9]  have  discussed  in  great  detail  a  non-Gaussian 
model.  On  the  other  hand,  R  ef  erenses  (6)  and  (10)  indicate 
that  filtered  wind  speeds  (particularly  with  a  band-pass  5-50 
Hz  filter)  is  reasonably  Gaussian  in  character.  In  the 
present  work,  the  Gaussian  model  is  used  exclusively. 

The  statistical  quantities  of  most  frequent  interest  are 
the  mean,  variance,  probability  density,  autocorrelation 
function,  power  spectral  density  and  crossing  frequency.  The 
power  spectral  shapes  usually  assumed  are  those  proposed  by 
von  Karman,  Kaimal  and  Teunissen.  The  crossing  frequency 
with  these  three  spectral  expressions  will  be  studied  in 
more  detail  later. 


form  characterized  by  three  parameters,  the  gust  vel¬ 
ocity  standard  deviation^,  the  scale  length  (L)  ,  and 
the  averaged  speed  (0)  . 

3)  Each  of  the  three  gusts  components  is  a  Gaussian 
process. 

The  assumptions  of  the  Gaussian  model  make  it  possible 
to  calculate  the  statistics  of  any  structure  response  for 
each  region  of  turbulence  as  functions  of  the  parameters  L 
and  o' of  that  region.  Thus,  the  assumptions  of  structure 
linearity  and  the  Gaussian  nature  of  turbulence  permit  their 
evaluation  with  a  minimum  of  difficulty.  These  results  will 
be  dependent  upon  the  assumed  values  of  L  and  <f  along  with 
the  characteristics  of  the  structure  dynamics. 

The  continuous  model  is  used  primarily  to  evaluate 
response  statistics  for  selected  conditions  in  continuous 
turbulence.  For  this  application  both  the  scale  length  and 
standard  deviation  of  the  turbulence  as  well  as  the  struc¬ 
ture  characteristics  are  fixed  at  values  representative  of 
certain  conditions  which  could  reasonably  be  expected  to 
occur  in  service. 

The  three  assumptions  made  above  need  some  reconsidera¬ 
tions.  The  first  assumption  stated  that  the  turbulence  is 
stationary  and  homogeneous  and  the  dimensions  of  structures 
are  much  smaller  than  the  scale  lengths  of  the  gusts  compo¬ 
nents.  These  conditions  are  not  always  satisfied.  The  second 
assumption  of  the  power  spectral  model  is  valid  at  low  alti¬ 
tudes,  but  is  not  so  certain  at  high  altitudes,  the  problem 
made  being  more  difficult  by  a  lack  of  data  at  the  latter 
altitudes.  The  third  assumption  stated  that  the  turbulence 
is  a  Gaussian  process.  However,  some  experimental  data  indi¬ 
cate  that  the  Gaussian  turbulence  model  significantly  under¬ 
estimates  the  frequency  of  occurence  of  high  gust  velocities 
[Bef.  7].  Furthermore,  peak  gust  velocity  cumulative  prob¬ 
abilities  of  exceedance  observed  in  atmospheric  turbulence 


sufficient  damping,  that  each  encounter  with  one  of 
these  gusts  results  in  a  single  significant  response 
peaks. 

This  model  is  typically  sed  to  estimate  ultimate 
strength  requirements,  however,  it  is  not  a  satisfactory  one 
for  all  aspects  of  the  design  problem.  For  example,  its 
application  in  areas  such  as  structural  fatigue,  control 
system  performance,  or  even  extreme  responses  which  involve 
lightly  damped  modes  is  highly  guestionable.  In  most 
instances  turbulence  does  not  occur  as  discrete  gusts  but  as 
a  continuous  random  disturbance.  Furthermore,  this  model 
neglects  most  of  the  dynamic  characteristics  of  the  vehicle. 
Thus,  although  the  necessity  of  evaluating  the  responses  of 
a  proposed  structure  to  discrete  gusts  is  still  recognized 
as  an  important  part  of  the  design  procedure,  the  use  of  the 
discrete  model  to  calculate  most  response  statistics  has 
largely  given  way  in  recent  years  to  the  use  of  turbulence 
models  which  attempt  to  take  both  the  dynamic  characteris¬ 
tics  of  the  structure  and  the  continuous  nature  of  turbu¬ 
lence  into  consideration.  This  is  known  as  the  continuous 
or  power  spectral  model. 


B.  CONTINUOUS  HODEL 

The  basic  idea  of  continuous  model  is  that  atmospheric 
turbulence  can  be  represented  by  a  continuous  stochastic 
process  which  acts  as  a  disturbing  influence  on  the  struc¬ 
ture.  The  principal  assumptions  for  this  model  are 

1)  Each  encounter  of  an  structure  with  continuous  atmos¬ 
pheric  turbulence  can  be  modeled  as  a  deterministic 
linear  system  perturbed  by  three  independent  stati¬ 
onary  stochastic  processes,  which  represent  the  long¬ 
itudinal,  lateral  and  vertical  gusts  components. 

2)  The  spectral  density  of  each  random  process  has  a 
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Figure  2.1  Definitions  of  GUSTO  and  GUST1. 

The  discrete  model  treats  turbulence  as  a  series  of 
isolated  gusts  and  the  basic  assumptions  are 

1)  Atmospheric  turbulence  can  be  modeled  as  a  collection 
of  isolated  gusts  randomly  distributed  along  the 
structure. 

2)  Gusts  have  random  magnitude  but  fired  shape. 

3)  The  structure  is  a  deterministic  linear  system  with 
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EFFECTS  OF  CUT-OFF  FREQ. 


B.  CALCULATION  OF  C BOSSING  FREQUENCIES 


In  this  section  the  frequency  of  zero  crossings  are 
estimated  using  von  Karman,  Kaimal  and  Teunissen  spectral 
density  functions  with  various  filter  functions  which  we 
have  discussed  above. 

1.  Using  the  von  Karman  Spectral  Density  Function 

Introducing  the  von  Karman  spectral  density  function 
for  the  x-direction  Su(n)  into  Equation  (4. 2)  yields 


N 


ox 


[  Pn*  F  (n)dn 

J(1+70.8n3)5/* 

'  He  i 


J\c» 


F(n)  dn 


BHJ 1* 


For  numerical  calculation,  set  y=n/nc,  and  D=nci  /10nC| 
Equation  (4. 13)  transformed  to 


(4.13) 


and 


^  ~  T  rloP 

Nox-nci'  ' 


(4. 14) 


where  N„*  is  the  dimensionless  frequency  of  positive- si ope 
level  crossings  for  the  x  direction  NOJ(Lu/U.  Equation  (4.  14) 
is  for  the  ideal  band-pass  filter  function  of  D  decades 
filter  width. 

The  upper  integral  of  Eg.  (4.13)  will  not  converge  as 
nCa  goes  to  infinity  because  the  empirical  expression  for 
the  turbulence  spectrum  employed  does  not  account  for  the 
high-frequency  range  where  viscous  dissipation  damps  out 
turbulence  fluctuations.  One  can  achieve  convergence  by  use 
of  a  suitable  filter  function.  Introducing  Equation  (4. 9) 
into  Equation (4. 13)  yields 


N0* 


nc 

r*  y3  <3y 

/  f  dY 

^(H70.8ncV)5/<5(ny*)/ 

1  +  70. 8nf  y2  f*(1  +  y2) 

O  v 

(4.  15) 
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If  we  take  F2 (n) ,  here  called  the  sine  filter  function, 
instead  of  quadratic  filter,  with  the  cut-off  frequency 
n,  =0.44/1  which  Greenway  [Ref.  12]  has  used,  and  substitute 
into  Equation (4. 6) ,  we  find 

F2  (n)  =sinz  (0.  44*n/nc  )  /(0  .44  tc,  n/nc  )Z.  (4.  16) 

Introducing  Equation  (4. 16)  into  Equation  (4. 13)  and  setting 
y=0.44un/nc  yields 


f  ft* 

f.  sina(y)dy _ 

N0*  *  _££_  J  %fl  +  70.  8  (Dc  /0.44it)2ya] 
0.44TC  |  r00  [sin2  (y) /y2]  dy 

J[U  " 


(4.  17) 


f[l  +  70.8(nc  /0 . 44TC)2  ya] 


If  we  use  Fj  (n)  ,  which  is  the  Hanning  frequency 
window,  and  equate  cut-off  frequency  nc=0.375/T  as  derived 
by  Vinnichenko  et  al.  [Ref.  11],  then  substitution  into 
Equation  (4.  8)  gives 

F^(n)=  /  sin  (0. 375icn)\  Tt* _  (4.18) 

l  0.375iui  /[it.2-  (0.375icn)2]S' 


Substituting  Equation  (4. 1 8)  intc  Equation  (4.  13)  and  setting 
y=0.375nn/nc  and  =‘nc/0.37  51l  yields 


% 


r  sin2 (y) dy 


A 


nc*-y2)a 


£ 


[sin2  (y)  /y23  dy 


7(1+70.8  7/  y 3  f/6  (IT2 - ya )2 


(4.  19) 
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T&BLE  I 

Crossing  Frequencies (von  Karman  z- direct ion) 


£c 

IDBP 

QUAD 

SINE 

HANNING 

0.  100 

0.3972 

0.  1215 

0.0865 

17.4480 

0.  125 

0.4838 

0.  1116 

0.0955 

21 . 8074 

0.  150 

0.5709 

0.  1250 

0.1117 

26. 1670 

0.  175 

0.6585 

0.  1381 

0.  1232 

30.5267 

0.220 

0.8172 

0.  1604 

0.  1424 

38.3741 

0.  500 

1.8193 

0.2737 

0.2391 

87.2049 

1.000 

3. 6235 

0. 42  9  8 

0.3699 

174.4062 

2.000 

7.2394 

0.6765 

0.5744 

348. 8093 

5.000 

18.093  1 

1.  2365 

1.0383 

872.0234 

10.000 

36.1  84  7 

1.9555 

1.6267 

1744.0481 

15.000 

54.2766 

2.  5585 

2. 1230 

2616.0669 

22.000 

79.6055 

3.2989 

2.7307 

3836.9021 

50.000 

180.9212 

5.  6928 

4.6950 

8720.2461 

75.000 

271.3816 

7.4555 

6.1414 

— 

100.000 

36  1.8425 

9. 0288 

7.4326 

— 

150.000 

54  2.76  37 

11.8268 

9.7292 

—  — 

175.000 

633.  2241 

13. 1054 

10.7787 

— 

200.000 

723.6853 

14.3243 

11.7794 

— 

225.000 

814.  1465 

11.4933 

12.9392 

— —  — 

250.000 

904.6067 

16.6196 

13.6641 

— 

nc=cut-off  frequency 
IDBP=ideal  band-pass  filter 
QUAD= "quadratic-type"  filter 
SINE= "sine-type"  filter 
BANNING”  HANNING  filter 
IDLPs=ideal  low-pass  filter 


The  results  of  numerical  calculations  of  Egs.(4.14),  (4.15), 
(4.17)  and  (4.19)  are  given  in  table  I.  The  plots  for  above 
equations  are  shown  in  Figure  4.4.  The  plots  showing  how 
well  these  crossing  frequencies  co-incide  depend  on  the 
length  scale,  mean  wind  speed  and  filter  functions. 

The  von  Karman  spectral  density  function  for  the 
z-direction  is  given  by  Equation (3. 3) .  The  equation  of  the 
crossing  frequency  for  the  ideal  band-pass  filter  function 
can  be  obtained  by  substituting  Equation  (3. 3)  into 
Equation  (4. 2)  and  is  given  by 
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VON  KARMAN(X-DIRECTION) 


LEGEND 
a  VKQUADX 


+  VKIDBPX 
X  VKHANNINGJT 


Using  the  same  procedure  as  for  the  x-direction,  the 
z-direction  crossing  frequency  for  the  quadratic  filter 
function  is  found  to  be 


'  J_ 

P**(  1+755.  2nexya )  y2dy 

9  (1+283. 2ncV  f6  P+y2) 

A 

^(1+755.  2ncV)  dy 

( 1+283.25* y2  V*  (i  +y4) 

(4.21) 


and,  for  the  sine  filter  function,  if  one  defines 
nc /0. 44lt  =»fc  the  crossing  frequency  for  the  z-direction 
becomes 


N 


OH 


PCl+755.2  (n,i  yflsin*  (Y) 

(  1  +  283. 25cV)W/* _ 

PCI +755.  2  (ncy)3]Qina  (y)  /Y*jdy 
,J0  (1  +  283. 2»{c2y2)"/6 


(4.  22) 


Introducing  filter  function  F3  (n)  into  Equation  (4. 20)  and 
setting  y=0.375itn/nc  and  nc /0.  375n.=^c  yields 


TABLE  II 

Crossing  Frequencies (von  Karman  z-direction) 


A* 

He 

IDBP 

QUAD 

SINE 

HANNING 

0.  100 

0.3812 

0.  121 

0.0816 

17. 4448 

0.  125 

0.4686 

0.  136 

0.0934 

21.8044 

0.  150 

0.5567 

0.  150 

0. 1044 

26. 1641 

0.  175 

0.6454 

0.163 

0.1144 

30.5238 

0.220 

0.8059 

0.  185 

0. 1320 

38.3716 

0.500 

1.8137 

0.297 

0.2188 

87.2036 

1.000 

3.6207 

0.452 

0.3375 

174.4053 

2.000 

7.2380 

0.696 

0.5238 

348.8093 

5.000 

18.0925 

1.248 

0.9458 

872.0229 

10.000 

36.  1344 

1.951 

1.4870 

1744.0461 

15.000 

54.2764 

2.  535 

1. 9407 

2616.0693 

22.000 

79.6054 

3.  246 

2.4978 

3836.8989 

50.000 

180.9214 

5.493 

4.2994 

8720.2266 

75.000 

271.  3823 

7.  127 

5.6256 

—  — 

100.000 

361.  84  30 

8.627 

6. 8094 

— 

150.000 

542.7634 

11.129 

8.9150 

— 

175.000 

633.2231 

12.330 

9.8773 

—  — 

200.000 

723.6838 

13.474 

10.7948 

— 

225.000 

814.  1440 

14.572 

11.6748 

— 

250.000 

904.6057 

15.63C 

12.5228 

— 

i - j 

The  results  of  numerical  evaluations  of  Egs.  (4.20)  , 
(4.21)  ,(4.22)  and  (4.23)  are  listed  in  the  Table  II.  The 
plots  for  above  equations  are  shown  in  Figure  4.5. 

For  comparison,  the  crossing  frequencies  with  von 
Karman  spectral  expression  for  the  x  and  z  directions  are 
plotted  in  Figure  4.6.  The  crossing  frequency  for  the  z 
direction  is  less  than  the  x  direction  except  ideal  band¬ 
pass  filter  case,  in  which  the  frequencies  are  identical  for 
both  directions. 

The  crossing  frequency  plots  presented  in  Figures 
4.4  -  4.6  offer  a  detailed  view  of  how  the  crossing 
frequency  depends  on  the  filter  function.  Thus,  in  designing 
a  system,  the  type  of  the  filter  has  to  be  chosen  carefully 
to  match  the  behavior  of  the  particular  system. 
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VON  K A R M A N (Z - DI R ECT ION) 


io  io-  io“  id 

CUT-OFF  FREQUENCV(nc) 

Figure  4.5  Crossing  Frequencies  with  von  K 
Spectral  Expression  for  the  z-direction 


arman 


fisifia  the  Ka imal  Spectral  Density  Function 

From  Eguation  (3.4)  ,  the  asymptotic  spectral  behavior 
can  be  written  as 


nS«  (n)  q  f  (f/fe)V3  ;  f>>f0 

#2*  l  0 . 164  (f/f0  )  ;  f<<f0 


(4.24a) 

(4.24b) 


VON  KARMAN(X  AND  Z  DIRECTION) 


Figure  4.6  Crossing  Frequencies  with  von  Karman 
Spectral  Expression  for  the  x  and  z -direction. 

The  expression  for  the  low  frequency  behavior  of  the 
spectrum  is  given 

•°S‘  ^  U/l)  f 

«I  (4.25) 

where  !<•  represents  the  length  scale.  Equating 

Equation  (4.25)  and  Equation  (4.  24b)  becomes 

=  0.  041z/(f.  )*.  (4.26) 
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.i. 


Por  the  x-direction  only 
L  n  =  0.04  1z/f0-x 


Setting  y=e*  in  Equation (4.  3 1)  gives 


No*  =nc 


e 3*  dx 


Vz 


[1+33.641  (nce*fej  (1+e*x) 

1  e*  dx _ 

[1+33.641  (ncex)%J(1  +  e25C) 


(4.32) 


If  we  introduce  the  sine  filter  function  into 
Equation  (4.29)  ,  the  crossing  frequency  becomes 


N 


ox 


sin*  (0.44TCD/nc)n2dn 
(1+33.641n*/3)  (0.44iC/ncf 
00  sin*  (0. 44nn/nc  )  dn 
(1+33.641nff/3)  (0.4 4it/nc  f 


(4.33) 


Substituting  0.44jtn/nc  for  y  in  Equation  (4. 33)  ,  and  define 
nc/0. 4411=^  /  yields 


N 


ox 


2 

sin  (y)  dy 
1+33.641  (yfcY)5/3 
"[sin2  (y)  /ya]  dy 
1  +  33.641  (ncy)5/3 


(4.34) 


For  the  Hanning  filter, 
Equation  (4. 28)  ,  setting 
nc/0.375tl=^c  ,  yields 


substituting  Equation  (4 . 1 8)  into 
y=0 .3  75]tn/nc  and  define 


i  , 

.  i 


1 
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f  r. _ 

N0^=ncJ  Jo  [1*33.641  (ft 


(y)£y 


y)f/3]  (i^-y2)* 


fsin2 (y) /y2] dy 


(1+33.  641  (ftcy)%]  (Tt2  -y2)2 


(4.35) 


TABLE  HI 

Crossing  Frequencies (Raimal  x-direction) 


Sc 

IDBP 

QUAD 

SINE 

HANNING 

0.  100 

0.4056 

0. 134 

0.0868 

17.4532 

0.  125 

0.4938 

0.  152 

0. 1001 

21.8126 

0.  150 

0.5819 

0.  169 

0.1126 

26.  1  720 

0.  175 

0.6702 

0.  184 

0. 1243 

30.5316 

0.220 

0. 8295 

0.210 

0.1441 

38. 3790 

0.500 

1.8305 

0.340 

0.2445 

87.2088 

1.000 

3. 6321 

0.518 

0.3811 

174.4094 

2.000 

7.2456 

0.  799 

0.5943 

348. 8120 

5.000 

18.0968 

1.434 

1.0744 

87  2.0224 

10.000 

36. 1872 

2.  249 

1.6879 

1744.0486 

15.000 

54.2786 

2.931 

2.1987 

26  16.0667 

22.000 

79.6070 

3.  769 

2.8249 

3836.9033 

50.000 

18 0.9221 

6.  471 

4.8291 

8720.2461 

75.000 

271.3823 

8.464 

6.2923 

— 

100.000 

361.  8423 

10.243 

7.5860 

— - 

150.000 

54  2.  76  34 

13.324 

9.8619 

— 

175.000 

633.  2239 

14.733 

10.8962 

— - 

200.000 

723.6846 

16.102 

11. 9076 

- — 

225.000 

814. 1450 

17.414 

12.8333 

— 

250.000 

904.6052 

1  8.678 

13.7017 

— 

The  results  of  numerical  calculations  of  Eqs.  (4.29), 
(4.  32),  (4.34)  and  (4.35)  are  given  in  Table  III  and  the 
plots  for  above  equations  are  shown  in  Figure  4.7. 

It  is  clear  from  the  Eg.  (3.4)  and  Eg.  (4.26)  that, 
for  the  Kaiaal  spectral  density  function,  the  difference  in 
crossing  frequencies  between  the  x  and  z  directions  is  only 
in  the  length  scale  .  Thus  the  plots  for  the  z-direction 
are  exactly  same  as  the  x-direction. 
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Figure  4.7  Crossing  Freguencies  with  Kaimal 
Spectral  Expression  fcr  the  x-direction. 

3.  Osing  the  Teunissen  Spectral  Density  Function 

The  modified  Kaimal  spectral  density  expressions, 
which  were  derived  by  Teunissen,  are  given  in  Egs.  (3.5), 
(3.6)  and  (3.7).  These  eguations  were  derived  by  adding  a 
scaling  parameter  to  each  of  the  Kaimal  spectral  eguations 
and  evaluating  these  parameters  from  the  observed  data. 
These  eguations  are  function  of  friction  velocity  u^  .  In 
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order  to  estimate  the  crossing  frequency  as  we  have  done 
earlier,  we  need  to  change  that  equations  as  a  function  of 
variance.  In  E.S.D.U.  [Ref.  5]#  the  relations  between  <3^ 


and  are  defined  as 

O.40£  /u*=  {<£  /U)  log  (z/2e)  (4-  36) 

Teunissen  derived  the  average  values  of  /u^,  given  by 

<ru/u*  =  2.  84  (4.37a) 

c^/u*  =  1.92  (4.37b) 

<5^/u*  =  1.27  (4.3  7c) 

From  Equation  (4 . 36)  and  Equation  (4.  37)  ,  one  can  get 

u*  (u-component)  =  0.124  (fu  (4.38a) 

u£  (v-component)  =  0.271  (4.38b) 

uj  (w-component)  =  0.6  20  (4.38c) 


and  substituting  Eg.  (4.38)  into  Eqs.(3.5),  (3.6)  and  (3.7) 

yields 

nSu  (n)  _  1 3n  (4.39) 

<Tua  (0.44  +  33n)5’/3 

nSy(n)  _  4. 6 In  (4.40) 

<TV2  ~  (0.38+9. 5n)S/3 

nSw(n)  1. 24n  (4.  41 ) 

o'*2  0.44  +  5.  3n*/s 

Using  methods  similar  to  the  above,  the  crossing  frequency 
for  the  x-direction,  with  the  ideal  band  pass  filter, 
becomes 


As  a  consequence  of  the  results  studied  herein,  it  is 
suggested  that,  in  future  work,  careful  consideration  should 
be  given  to  the  type  of  filtering  used  to  model  structural 
behavior  and  determine  more  accurately  the  fatigue  life  of  a 
system. 


?.  CONCLUSIONS 


This  study  has  estimated  the  crossing  frequency  based  or. 
the  three  spectral  density  functions,  which  assumes  turbu¬ 
lence  to  be  a  homogeneous,  stationary  Gaussian  process,  with 
five  different  filter  functions  for  the  x  and  z  directions. 

This  estimation  has  shown,  for  the  x  direction,  that  the 
Kaimal  spectral  function  gives  about  one  percent  higher 
crossing  frequency  than  the  von  Karman,  and  forty  percent 
higher  than  the  Teunissen  one,  at  2^  =  100  for  quadratic-type 
filtering;  for  the  ideal  low-pass  and  sine-type  filters,  the 
difference  between  the  Kaimal  and  von  Karman  spectral  func¬ 
tion  is  only  a  few  percent,  tut  for  the  Teunissen  one  is 
more  than  forty  percent  than  for  the  other  spectrums. 

For  the  z  direction,  the  Teunissen  spectral  function 
gives  about  ten  percent  higher  crossing  frequency  than  the 
Kaimal  one  and  about  twenty  percent  higher  than  the  von 
Karman  one,  at  nc=100  for  the  quadratic- type,  sine-type  and 
low-pass  filtering.  However,  this  difference  becomes  to  zero 
as  non-dimensional  cut-off  frequency  goes  to  zero. 

The  crossing  frequencies  for  band-pass  and  Hanning 
filtering  are  only  a  fraction  of  a  percent  different  for  the 
three  spectral  functions. 

The  present  study  shows  that  Hanning  and  band-pass 
filtering  predicts  higher  crossing  frequencies  than  low-pass 
filtering.  This  difference  is  greatly  amplified  with  the 
higher  frequency. 

If  some  site  statistics  are  available,  as  for  the 
example  calculation,  then  the  number  of  crossings  that  the 
wind-speed  fluctuations  exceed  their  zero  value  per  a 
certain  period  can  be  computed  and  an  estimation  of  fatigue 
life  can  be  obtained. 
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Figure  4.13  Values  of  the  Surface  Roughne 
(Reproduced  fros  ESDD  74031). 


hness  Parameter  z 


. 
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TABLE  VI 

Example  Calculations  of  the  Crossing  Frequency 


REMARK 

Nox 

Nox 


REMARK 

ifoz 

Noz 


REMARK 


Nox 
ox 
oz 
Noz 


if 


REMARK 

Nox 

Nox 

NOZ 

NOZ 


von  Karmantx 

:-d  ires 

tion) 

I  EBP 

ID  IP 

SINE 

75.0 

10.2 

14.6 

12.0 

17.857 

2.  428 

3.476 

2.  857 

(2cu=210.0 

,  0/Lu=0 .238) 

von  Karman ( 2 

;-direction) 

I  EBP 

ID  IP 

QUAD 

SINE 

4.70 

1.  50 

2.40 

1.  80 

18.076 

5. 769 

9.231 

6.923 

(8ew  =  1  » 

U/Lv=3 

.846) 

Kaimal  (x  and 

.  z-direction) 

I  EBP 

I  DIP 

QUAD 

SINE 

75.0 

10.9 

16.50 

12.2 

17.857 

2.  595 

3.928 

2.905 

4.70 

1. 70 

2.65 

2.05 

18.076 

6.  538 

10.  192 

7.885 

Teunissen lx 

and  zz 

direction) 

I  EBP 

I  DIP 

QUAD 

SINE 

75.0 

6.  08 

9.30 

6.50 

17.857 

1. 448 

2.214 

1.548 

4.70 

2.  10 

3.30 

2.40 

18.076 

8. 076 

12.  692 

9.230 

zo=0.001.  The  length  scale  parameters  are  given  in  Figure 
4.12,  and  are  given  as  Ltt=105,  Iv=36  and  1^=6. 5. 

The  next  step  is  calculation  of  the  dimensionless 
cut-off  frequencies,  i.e.. 


ncu  = 

nLu/D 

=  210 

(4.  50) 

Scv  = 

nLv/0 

=  72 

*cw  = 

nL^/tJ 

=  13 

The  crossing  frequencies  are  calculated  by  the  use  of 
the  given  crossing  frequency  vs  cut-off  frequency  graphs, 
with  dimensionless  cut-off  frequencies  which  are  given  in 
Eg.  (4.50). 

The  final  step  is  to  convert  the  dimensionless  crossing 
frequency  to  a  dimensional  frequency  by  multiplying  the 
factor  U/L^,  i.e.. 

No*  =  Nox  0/Iu  (4.51) 
Nol  =  No*  O/ljy 

The  results  of  calculation  cf  the  crossing  frequency  for 
a  given  turbulence  condition  are  given  in  Table  71. 


CROSSING  FREQUENCY(No*) 


© 

CO 


VON  KARMAN,  KAIMAL&  TEUNISSEN(Z) 


Figure  4.11  Comparison  of  Ton  Karman.  Kaimal  and  Teunissen 
Spectral  Density  Function  for  the  z-direction. 


Suppose  we  have  a  site  for  which 
z  =  20m 
nc=  50Hz 

0  =  25m/sec  and  the  turbulence  condition  is  that 
for  an  off-sea  wind  in  a  coastal  area. 

With  this  choice  of  turbulence  condition,  the  surface 
roughness  parameter  ze  can  be  found  in  Figure  4.13,  and  is 


VON  KARMAN,  Kf IMAL&  TEUNISSEN(X) 


Figure  4.10  Comparison  of  von  Karman.  Kainal  and  Teunissen 
Spectral  Density  Function  for  the  x-direction. 

The  dimensionless  crossing  frequency,  Ne^ ,  is  a  function 
of  nc  (*nL,i  /U) .  The  length  scale  parameter  and  surface 
roughness  parameter  z0  are  giver  by  E.S.D.O.  [Ref.  5].  For 
convenience  these  data  are  included  here  as  Figure  4.12  for 
length  scale  parameters,  and  Figure  4.13  for  the  surface 
roughness  parameters. 
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CUT-OFF  FREQUENCY(ilc) 


Figure  4.9  Crossing  Freguencies  with  Teunissen 
Spectral  Density  Function  for  the  z-direction. 

C.  EXAMPLE  CALCULATIONS 

In  the  previous  section  of  this  study,  we  derived  equa¬ 
tions  for  the  crossing  frequency.  The  purpose  of  this 
section  is  to  illustrate  the  application  of  these  equations 
to  a  given  turbulence  condition,  and  to  estimate  the 
crossing  frequencies  for  a  low-pass  filter. 


TABLE  V 

Crossing  Freguencies (Teunissen  z-direction) 


£c 

I DBP 

fiUAD 

SINE 

IDLP 

0. 100 

0.4454 

0.178 

0. 1098 

0. 0559 

0.  125 

0.5367 

0.20  1 

0.1258 

0.0691 

0. 150 

0.6263 

0.222 

0.  140  8 

0.0818 

0.  175 

0.7151 

0.241 

0.1551 

0.0  94  1 

0.220 

0.8740 

0.274 

0.1794 

0. 1153 

0.  500 

1.8651 

0.436 

0.3047 

0. 2268 

1.000 

3.6560 

0.660 

0.4758 

0.3783 

2.000 

7.2610 

1.007 

0.  7414 

0.6104 

5.000 

18. 1053 

1.789 

1.3353 

1.1232 

10.000 

36. 1926 

2.791 

2.0917 

1.7725 

15.000 

54.2827 

3.629 

2.7242 

2.3144 

22.000 

79.6102 

4.658 

3.5000 

2.9785 

50.000 

180.9239 

7.973 

5. 9955 

5. 1216 

75.000 

271.3835 

10.410 

7.8197 

6.6982 

100.000 

361. 8437 

12.  591 

9.4420 

8. 1050 

150.000 

542.7644 

16.4  70 

12.3095 

10.6064 

175.000 

633.2249 

18. 242 

13.6158 

11.7493 

200.000 

723.6853 

19.9  32 

14.86  06 

12.8388 

225.000 

814.  1458 

21.553 

16.0474 

13.8837 

250.000 

904.6062 

23.  1 15 

17.1816 

14.8905 

For  comparison,  the  crossing  freguencies  with  the 
von  Karman,  Kaimal  and  Teunissen  spectral  expressions  for 
the  x-direction  with  various  filter  functions  are  plotted  in 
Figure  4.10.  As  shown  in  Figure  4.10,  the  crossing 
freguencies  predicted  by  the  Kaimal  spectral  expression  are 
higher  than  for  the  others.  In  Figure  4.11,  the  crossing 
frequencies,  with  all  three  spectral  expressions  and  the 
z-direction,  are  compared.  The  figure  characteristics  are 
similar  to  those  for  the  x-direction,  except  that  the 
crossing  freguencies  are  lower. 
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If  one  introduces  the  sine  filter  function,  the  equation  for 
the  crossing  frequency  becomes 


*og 


ft 

JoJ! 

r"  sin2  (0. 4 
J  (0.44+5. 3n 


sin3  (0.44nn/nc)  na  dn 
(0.44+5.  3n%)  (0. 44Un/nc  f 

44itn/nc)  dn _ 

*4)  (0.44nn/nc)zJ 


(4.47) 


r 

r°°  2 

/  sin  (y)dy 

or 

N0*  =nc  - 

l  0.44+5.3  incY)*/3 

u 

'"[sin2  (y)  /y2]  dy 
0-44+5.3  <ffc  y)5/^ 

l/z 


(4.48) 


j£  one  introduces  the  quadratic  filter  function,  the  egua~ 
tion  for  the  crossing  frequency  becomes 


M 


Y2) 


r  Yza  y _ 

Je[0.44+5.3  (ncy)^]  (1  + 

r* _ dy _ _ 

j  [0.44+5.  3  tfcyf/3]  (i  +  yz) 


A 


(4. 49) 


The  results  of  numerical  calculations  of  Egs.  (4.46) , 
(4.48)  and  (4.49)  are  given  in  Table  V.  The  plots  for  above 
equations  are  shown  in  Figure  4.9. 


CUT-OFF  FREQUENCY(nc) 


Figure  4.8  Crossing  Frequencies  using  Teunissen 
Spectral  Expression  for  the  x-direction. 


TABLE  I¥ 

Crossing  Frequencies  (Teunissen  x-direction) 


Sc 

I  DB£ 

quad 

SINE 

IDLP 

0. 100 

0.3722 

0.060 

0.0425 

0.0334 

0.  125 

0.4628 

0.069 

0.0494 

0.0394 

0.  150 

0.5534 

0.078 

0.0559 

0.0449 

0.  175 

0.6439 

0.086 

0.0620 

0.0502 

0.220 

0.8068 

0.099 

0. 0724 

0.0590 

0.  500 

1.8201 

0. 16S 

0.1251 

0.1044 

1.000 

3.6294 

0.266 

0.  1983 

0. 1672 

2.000 

7.2478 

0.420 

0.3142 

0.2664 

5.000 

18.  1031 

0.768 

0.5754 

0.4915 

10.000 

36. 1952 

1.213 

0.9083 

0.7801 

15.000 

54.  2873 

1.581 

1.  1846 

1.0219 

22.000 

79.6162 

2.0  35 

1.5228 

1.3188 

50.000 

180.9321 

3.521 

2.  6002 

2.2786 

75.000 

271.3931 

4.612 

3.3801 

2.9851 

100.000 

361.8535 

5.586 

4.0689 

3.6157 

150.000 

542.7749 

7.266 

5.2816 

4.7372 

175.000 

633.2354 

8.053 

5.8554 

5.2496 

200.000 

723.6956 

8.804 

6. 2870 

5.7380 

225.000 

814.1565 

9.525 

6.8032 

6.2065 

250.000 

904.6  169 

10.219 

7.301  1 

6.6579 

The  results  of  numerical  calculations  of  Egs. (4.42)  , 
(4.44)  and  (4.45)  are  given  in  Table  IV  and  the  plots  for 
these  equations  are  shown  in  Figure  4.8. 

The  modified  Kaimal  spectral  expression  for  the 
z-direction  is  given  by  Eg.  (4.41).  If  we  put  ideal  band  pass 
filter  function  and  spectral  expression  for  the  z-directon 
into  Eq.  (4.2)  it  yields 


N 


oz 


r  n  dn  / 

/r"‘2a  n 

J  0.44+5. 3n5/1»  / 

J  0.44+5.3n5^ 

or 


^ox  = 

r"ta  nxdn 

r"‘2  dn 

J  (0.44*33nfV 

L  nci  ' 

J  (0.44+33n)5/3 

Tlei  ^ 

»o*  * 

nc. 

r°cl  y2dy 

/r"61  dy 

Vi 


nci 


nci 


(4.42) 


If  one  uses  the  sine  filter  function,  with  cut-off 
frequency  nc=0.44/T,  the  crossing  frequency  becoaes 


N 


ox 


sin*  (0. 44Tin/cc  )  n*  dn 
(0.  44*33n)5/*  (0.44nn/nc  )2 
88  sin*  (Q.  44itn/nfc )  dn 
(0.44+33n)5/*  (0 .44trn/nc  )* 


(4.43) 


Setting  y=0.44itn/nc  into  Eg.  (4.43),  and  define  nc /0 . 4 4Tt=  i|c  , 
yields 


* 

► 


sin2  (y) dy 
(0. 44*33  r[cy)% 


/  fCsin2  (y)/y3]  dy  Y1 

/  ]q  (0.44*33^cy)H  f  (4.44) 


For  the  quadratic  filter  function,  becomes 


09 

h 


yady 


(0.44*3 3 nc  Y)f/i  ( 1»Y*) 


f  - 

J  (0. 


dy 


(0.44*33nc  y)f/s  (1*ya) 


Vi 


(4.45) 
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